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consideration  In  the  selection  of  specific  systems  for  study,  particularly 
multilayer  Integrated  circuit  applications. 

In  this  period,  experiments  with  the  new  cryopumped  system  were  performed. 
Reactive  sputtering  research  focused  on  a  detailed  study  of  S10x  and  a  more 
broad  guage  study  of  M0X  films  for  a  number  of  metals  M  (M  =  Si,  W,  At,  Zn,  Zr, 
Mg)  on  3H  diameter  wafers.  High  deposition  rates,  a  wide  range  of  stoichiometry 
and  either  amorphous  or  crystalline  films  were  formed.  The  cleanliness  of  the 
new  system  Is  orders  of  magnitude  better  than  the  old  system.  Multiple  source 
deposition  studies  focused  on  $1Cx  formation.  X-ray  analysis,  refractive  Index 
and  bandgap  data  all  confirm  the  formation  of  SIC  films.  Good  control  of  the 
sputter  deposition  of  these  films  has  been  demonstrated. 

In  the  theoretical  area,  the  constraints  placed  on  closed  n-rlng  assemblies 
(S1n03n)  In  S102  are  being  assessed  to  Identify  those  constraints  which  dlfferen 
tlate  the  structure  between  the  crystalline  and  vitreous  forms.  Interactions 
between  various  interstitial  species  and  dimer  (S1207)  units  or  tetramer 
(SI4O1?)  nn9s  ™  s^lica  have  been  carried  out  using  CNDO/2.  Great  progress 
has  been  made  In  evaluating  the  total  energy,  Ey,  at  0°K  for  the  following 
Interstitial  species.  Si,  Si2+,  0,  0*.  02+,  02>  H20,  OH"  and  H+. 
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INTRODUCTION 


During  this  third  six-month  period  of  the  program,  the  experimental 
emphasis  focused  on  reactive  sputtering  of  a  variety  of  oxides  with  continued 
in-depth  study  on  SiOx  plus  continued  work  on  the  two-source  deposition  of 
SiCx  films.  In  all  areas,  meaningful  progress  has  been  made  and  it  is 
especially  significant  to  see  how  well  these  techniques  can  be  extended  to 
other  oxide  materials.  The  theoretical  emphasis  continued  on  both  the 
structural  and  the  energetic  aspects  of  the  SiOx  material  and  significant 

gains  are  reported  in  those  sections  of  this  report. 

A.  Experimental  Results 

1.  Reactive  Sputtering 

(a)  SiOv 

In  the  previous  report^1),  some  experimental  evidence  indicated 
that  oxygen  introduced  into  the  chamber  at  the  substrate  surface  was  back- 
diffusing  to  the  silicon  sputtering  source  and  contaminating  the  surface.  In 
addition,  CV  plot  data  suggested  that  hydrocarbon  contamination  was  present  in 
the  Si 0X  films. 

During  the  present  period,  construction  was  completed  on  a  sputtering 
system  designed  to  deposit  uniform  films  on  3"  diameter  substrates  (see  Fig. 
la).  This  system  provides  several  advantages  over  the  one  described  in 
previous  reports.  Specifically: 

(1)  A  helium  cryopump  has  eliminated  the  problem  of  l\ydrocarbon 
contamination  which  resulted  from  back-diffusion  of  oil  in  the  previous 
system. 

(2)  The  Increased  pumping  speed  of  this  system  allows  for  control led- 
stoichiometry  syntheses  at  a  much  lower  reactive  gas/argon  ratio,  which 


reduces  the  previous  tendency  toward  target  oxidation  at  high  Og 
pressure. 

(3)  An  oxygen  dispersion  ring  (see  Fig.  lb)  further  isolated  the  target  from 
the  reactive  gas  while  providing  a  uniform,  radially  symmetric  flow  of 
oxygen  to  the  substrate  surface. 

(4)  The  use  of  a  3"  target  has  increased  the  area  of  film  uniformity.  The 
system  illustrated  in  Fig.  la  will  support  a  target  up  to  6"  diameter, 
should  this  become  necessary. 


Electron  microprobe  analysis  has  confirmed  our  anticipatd  correlation  of 
refractive  index  with  stoichiometry  by  coinciding  with  previously  measured 
SIPOS  films  (see  Fig.  2).  this  implies  that  low  temperature  fabrication 
produces  films  with  a  microstructural  integrity  at  least  equal  to  that  found 
with  the  high  temperature  CVD  process. 

The  present  emphasis  is  on  determining  the  optimum  combination  of 
sputtering  parameters  (e.g..  Si  deposition  rate,  oxygen  pressure,  argon 
pressure)  which  will  produce  maximum  film  uniformity.  Since  uniformity  is 
affected  by  the  nature  of  the  oxygen  flow  over  the  substrate,  it  is  important 
to  have  a  wide  effective  range  of  oxygen  pressure  available  for  a  particular 
synthesis.  Fortunately,  there  is  a  well  behaved  relationship  between  gas 
pressure,  sputtering  rate  and  resulting  oxide  stoichiometry.  Figure  3 
illustrates  a  typical  experiment  in  which  oxygen  pressure  and  RF  power  are 
varied  independently  to  cover  a  full  stoichiometric  range.  At  500  W,  we  note 


that  the  film  stoichiometry  varies  linearly  with 
stoichiometry  can  be  simply  changed  by  changing 


so  that  the  film 
Likewise,  at  constant 


Pn  ,  the  film  stoichiometry  can  be  simply  changed  by  changing  the  sputter 
u2 


power. 
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The  radial  symmetry  of  our  present  samples  allows  film  uniformity  to  be 


easily  checked  using  ellipsometric  techniques  and  continued  improvement  is 
being  made  in  this  area.  To  illustrate  the  film  uniformity  in  the  new  system. 
Fig.  4  indicates  results  obtained  on  an  SiOj^  film.  We  note  that  the  index 
of  refraction,  ,  is  quite  uniform  out  to  a  1"  radius  from  the  center  of 

the  wafer  and  then  falls  off  seriously  at  ~  0.25"  from  the  outer  periphery  of 
the  wafer.  Over  this  same  film  distance,  the  wafer  thickness  is  found  to  be 
constant  within  ~  15%  .  The  next  step,  here,  will  be  to  optimize  the 
position  of  the  oxygen  flow  relative  to  the  substrate  surface,  thus  allowing  a 
much  more  uniform  oxygen  incidence  rate. 

Procedures  have  now  been  established  for  routinely  measuring  electrical 
properties  of  each  of  our  samples.  Monitoring  the  resistance  and  capacitance- 
diode  characteristics  will  ensure  that  the  oxide  layers  are  passivating  and 
that  the  interface  integrity  is  uniform.  These  measurements  will  lead 
directly  to  the  more  thorough  evaluation  of  film  electrical  behavior  (A.C.) 
planned  for  the  future. 

(b)  Extension  to  Other  Oxide  Systems 

During  this  period,  a  series  of  exploratory  experiments  into  the 
synthesis  of  the  oxides  of  aluminum,  zirconium,  magnesium,  and  zinc  were 
performed,  the  purpose  of  these  experiments  was  to  determine:  (1)  Can  this 
process  be  applied  to  other  oxide  systems?  (2)  What  characteristics  in  the 
synthesis  of  these  other  oxides  differ  from  those  of  SiOx.  (3)  Do  the  results 
suggest  that  this  process  can  be  applied  to  the  synthesis  of  other  compounds; 
i.e.,  to  those  containing  N,  C,  Cl  ,  F,  etc?  (4)  What  are  the  limitations  of 
this  process?  In  the  following,  the  results  of  these  experiments  will  be 
outlined  and  a  short  summary  included. 


The  results  of  these  experiments  are  summarized  in  Table  I.  It  was  found 
possible  to  synthesize  oxides  of  all  materials  attempted  except  magnesium, 
where  difficulties  with  the  sputtering  target  limited  the  experiments.  Only 
two  of  the  materials  deposited  were  found  to  form  crystalline  films  under  the 
conditions  explored  (zinc  and  zirconium).  ZnO  was  formed  with  control  of 
stoichiometry  and  conductivity.  Zr02  was  stabilized  in  the  high  temperature 
cubic  form,  with  an  apparently  amorphous  phase  appearing  as  the  oxide  was 
reduced  by  synthesis  at  lower  oxygen  pressures.  All  other  films  were 
amorphous  as  formed. 

Also  included  in  this  table  are  the  heats  of  decomposition,  AH,  of  these 
materials  at  ambient  conditions  and  $,  the  ratio  of  the  molecular  volumes  of 
the  oxides  to  that  of  their  metallic  components  in  their  elemental  forms. 

A*0_x 

Experiments  on  synthesis  of  AtOx  are  summarized  in  Table  II.  Both  theAt 
deposition  rate  and  the  oxygen  pressure  were  systematically  varied.  It  was 
observed  that  all  films  formed  were  non-stoichiometric  and  amorphous  (x-ray 
diffraction  analysis).  Composition  was  determined  by  electron  microprobe 
analysis  for  samples  197  and  199.  the  refractive  index  of  199  was 
nf  *  2.0  from  standard  ellipsometric  observations.  This  is  somewhat  larger 
than  o  -  A*203  but  consistent  with  a  higher  aluminum  concentration.  Sample 
197  was  examined  by  electron  microprobe  analysis  on  a  radius  of  the  sample; 
these  results  are  shown  in  Fig.  5.  Note  that  near  stoichiometric  oxide  is 
present  near  the  edge  of  the  sample  where  the  aluminum  deposition  rate  is 
smallest.  Also,  a  transition  from  clear  to  translucent  brown  was  observed  for 
0/At  ratios  less  than  approximately  1.25.  This  may  be  associated  with  color 
centers  or  aluminum  clusters  resulting  in  Mie  scattering  and  absorption. 
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The  variation  in  oxygen  to  aluminum  ratio  with  position  on  a  radius  of 
the  sample  can  be  interpreted  as  resulting  from  the  spatial  variation  in  the 
incidence  rate  of  aluminum  atoms  with  position,  which  is  characteristic  of  the 
magnetron  sputter  deposition  source  used.  If  one  assumes  that  the  oxygen 
incidence  rate  is  constant  over  the  sample  surface  and  that  the  sticking 
coefficient  or  reaction  rate  is  uniform,  the  thickness  profile  of  the  aluminum 
can  be  unfolded  from  these  data.  This  assumption  requires  that  the  number  of 
oxygen  atoms  incorporated  into  the  film  is  constant  and  independent  of 
position  on  the  substrate.  The  thickness  of  aluminum,  therefore,  varies  in 
ratio  from  1.0  at  10  urn  from  the  sample  center  to  0.8  at  20  mm  and  0.7  at  30 
mm.  This  thickness  variation  is  in  general  agreement  with  the  expected 
deposition  profile  of  the  aluminum,  implying  that  the  oxygen  gas  distribution 
is  sufficiently  uniform  for  synthesis  of  uniform  composition  films  if  the 
aluminum  incidence  rate  is  constant. 

These  results  also  demonstrate  that  the  stoichiometry  of  the  deposited 
films  can  be  systematically  varied  over  a  broad  range  by  control  of  the  rates 
of  incidence  of  the  oxygen  and  aluminum,  the  stoichiometry  is  a  sensitive 
function  of  the  relative  rates  of  incidence  requiring  control  of  this 
parameter  if  uniform  films  are  to  be  deposited. 

Possible  application  of  this  process  may  center  on  the  use  of  these  films 
as  passivation  layers  on  II I -V  compounds,  as  well  as  unique  passivation  layers 
for  more  exotic  materials. 


Synthesis  process  parameters  for  Zr/0  films  are  tabulated  in  Table  III. 
All  films  were  deposited  onto  Si (100)  single  crystal  substrates.  Analysis  was 
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performed  using  standard  powder  x-ray  diffractometry,  electron  microprobe 
analysis  and  ellipsometry. 

The  zirconium  deposition  rate  was  maintained  constant  in  these 
experiments  while  the  oxygen  pressure  varied  by  a  factor  of  four.  In  all 
cases,  the  films  were  near  stoichiometric  ZrOg.  The  film  formed  at  the  lowest 
oxygen  pressure  showed  a  slight  deviation  from  stoichiometry.  All  films  were 
observed  to  be  fine-grain  polycrystalline,  having  the  cubic  structure 
(aQ  *  5.01  A,  FCC).  This  cubic  structure  is  not  expectd  at  room  temperature, 
the  equilibrium  form  being  monoclinic.  The  cubic  form  is  generally  only 
observed  at  temperatures  above  ~  2300UC  in  equilibrium  and  thus  represents  a 
true  metastable  state.  Samples  formed  at  lower  oxygen  pressures  also  showed 
indications  of  a  broad  shallow  diffraction  maxima  which  could  be  associated 
with  either  an  amorphous  component  or  the  multiplicity  of  diffraction  peaks 
associated  with  orthorhombic  structure  Zr02* 

The  observation  of  a  high  reactivity  of  oxygen  with  the  zirconium,  as 
evidenced  by  the  insensitivity  of  the  film  structure  to  oxygen  pressure,  is 
consistent  with  the  common  use  of  zirconium  and  titanium  as  oxygen  getters. 

The  crystallinity  of  the  deposition  product  and  the  cubic  structures  are 
unique  and  indicate  that  material  specific  to  oxygen  sensor  applications  can 
be  deposited  in  thin  film  form. 

ZnOv 

Synthesis  process  parameters  for  Zn/0  films  are  tabulated  in  Table  IV. 

All  films  were  deposited  on  SI (100)  single  crystal  substrates.  Analysis  was 
performed  using  x-ray  powder  diffractometry,  ellipsometry  and  simple 
conductivity  tests. 
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In  all  cases  the  films  were  observed  to  have  the  hexagonal  crystalline 
structure  expected  for  ZnO.  Preliminary  measurements  of  the  lattice  parameter 
as  a  function  of  position  showed  no  dependence  upon  location  on  the  substrate 
and  therefore  on  stoichiometry  of  the  film.  Film  electrical  resistance  varied 
with  radial  position  from  the  center  of  the  sample  to  the  edge,  as  shown  in 
Fig.  6.  The  electrical  resistivity  of  ZnOx  is  a  function  of  stoichiometry  as 
reflected  in  these  results.  Film  composition  is  yet  to  be  characterized  by 
electron  microprobe  analysis. 

These  results  demonstrate  that  the  modified  experimental  apparatus  allows 
synthesis  of  oxides  of  controllable  stoichiometry  of  elements  such  as  Zn,  a 
very  high  sputtering  rate  material  having  a  very  stable  oxide.  Possible 
applications  include  transparent  conducting  films  for  windows  and  for  solar 
cell  contacts.  It  is  important  to  note  that  these  films  were  deposited  onto 
room  temperature  substrates  which  were  isolated  from  the  sputtering  plasma. 
These  are  important  characteristics  for  many  possible  applications  of  this 
material . 

Summary 

The  experimental  results  presented  in  the  preceding  discussion 
demonstrate  the  general  applicability  of  the  synthesis  technique  used  for  the 
formation  of  Si 0X  films.  This  process  allows  the  synthesis  of  whole  new 
classes  of  materials  onto  electronic  device  materials  which  can  be  expected  to 
act  as  passivating  layers  or  active  components  in  the  device.  Formation  of 
oxide  film  onto  room  temperature  substrates  at  deposition  rates  of  5  A/sec 
(18,000  A/hr)  with  control  of  stoichiometry  and,  in  some  cases,  of  the 
crystalline  structure,  represents  a  unique  capability.  This  processing  method 
makes  possible  the  synthesis  of  materials  which  have  not  previously  been 
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accessible  to  the  micro  electronics  community  for  evaluation  as  device 
components.  In  addition,  it  can  be  used  to  synthesize  Si 0X  films  onto 
substrates  held  at  low  temperatures  so  that  VLSI  structures  may  be  fabricated 
with  minimum  thermal  history. 

2.  Synthesis  of  SiC^Films  by  Two-Source  Deposition 

There  is  a  high  level  of  interest  in  the  semiconductor  device 
possibilities  of  B-SiC  ,  the  cubic  form,  the  temptation  to  generate  SiC  by  a 
process  compatible  with  silicon  technology  has  directed  this  part  of  the 
program  toward  low  temperature  (<  1000UC)  deposition  of  SiC  by  sputtering. 

The  system  used  has  a  combination  of  state-of-the-art  features  for  sputter 
deposition:  magnetron  sources  (plasma  decoupled  from  the  deposition  surface), 
substrate  heater,  cyro-pumped  chamber,  and  single,  dual  or  triple-source 
arrangement,  parallel  or  at  various  angle  to  the  substrate/deposition 
surface.  The  goal  is  to  achieve  deposition  of  large  areas  of  single  crystal 
B-SiC  by  a  sputtering  technique  at  a  temperature  with  typical  silicon 
processing  temperatures. 

The  set  of  experiments  described  here  indicates  that  Si  can  be  deposited 
on  Si  or  Az^Og  by  co-sputtering  silicon  and  carbon  targets.  SiC  films  having 
a  preferred  orientation  were  formed  at  700UC  on  Az^Og  and  at  750°C  on  Si 
substrates.  In  addition  to  investigating  the  structure  of  the  films,  various 
SiC  film  properties  have  been  measured.  The  refractive  index,  band  gap,  and 
stoichiometry  were  evaluated  using  standard  techniques.  These  analyses  will 
be  described.  The  present  experiments  are  a  continuation  of  the  dual-source 
depositions  described  in  the  previous  semi-annual  report^.  Deposition 
reproducibility  has  been  established  and  the  various  analyses  of  the  SiC 
indicate  that  B-SiC  has  been  formed.  The  next  set  of  experiments  will 
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involve  a  sputtering  configuration  that  should  enhance  the  surface  mobility  of 
Si  and  C,  resulting  in  SiC  films  of  greater  crystalline  perfection.  This 
configuration  and  existing  supporting  data  will  be  described  in  the  section  on 
future  direction. 


Sputter  Deposition  of  SiC 

First  the  sputtering  configuration  and  deposition  parameters  will  be 
reviewed.  In  Fig.  7  a  schematic  is  given  of  the  arrangement  inside  the 
sputtering  chamber.  Two  sources  were  employed;  a  graphite  source  was  operated 
in  tne  DC  mode  and  a  single  crystal  silicon  source  was  operated  in  the  RF 
mode.  The  source-to-source  separation  was  16.5  cm.  The  substrate  could  be 
heated  up  to  850°C  by  means  of  a  heat-shielded  quartz  lamp  configured  so  that 
the  substrate  was  heated  from  the  back  surface.  The  source-to-substrate 
distance  was  10-15  cm.  The  combination  of  the  magnetron  sources  and  this 
large  source-to-substrate  distance  were  considered  sufficient  to  keep  the 
sputtering  plasma  from  influencing  the  reaction  at  the  substrate  surface. 

Before  initiating  the  plasma,  the  sputtering  chamber  was  pumped  down  to 
2  x  10"®  Torr.  A  plasma  was  maintained  with  4  microns  (4  x  10"^  Torr)  of 
Argon.  As  indicated  in  Fig.  7,  a  monotonical ly  varying  distribution  of 
sputtered  material  is  expected  from  each  target.  With  the  targets  facing 
directly  downwards,  their  surfaces  parallel  to  the  substrate  surface,  an 
overlap  of  the  Si  and  C  distributions  is  evident  so  that  the  entire  range  of 
Si Cx  is  deposited  on  the  substrate.  One  end  of  the  deposition  is  carbon-rich, 
the  other  end  is  silicon-rich,  and  stoichiometric  SiC  is  somewhere  in  between. 

Two  types  of  substrates  were  used.  Si  and  M^y  The  ^  was  device" 
quartz  (100)  20n-cm  p-type  Si  which  had  been  cleaned  by  the  industry-standard 
"RCA  clean";  this  leaves  a  reproducible  (amorphous)  oxide  film  on  the  silicon 


surface.  The  A^O-j  was  SOS-quality  sapphire,  polished  on  the  side  used  for 
deposition,  and  cleaned  the  same  way  as  the  silicon.  The  A*203  is  not 
considered  to  have  any  type  of  amorphous  film  on  its  surface  and  should  act  as 
a  single-crystal  deposition  surface. 

Analysis  of  SiC  Films 

As  explained  at  the  ena  of  the  sputter  deposition  section,  there  is  a 
calculable  distribution  of  sputtered  material  from  each  target.  Thus,  the 
deposition  thickness  as  a  function  of  location  on  the  substrate  can  be 
anticipated  given  the  experimental  sputtering  rates  of  the  targets. 
Calculations  were  made  for  the  SiC  depositions  and  an  example  was  given  in  the 
last  semi-annual  report  (see  pages  10-11,  Table  IV,  and  Fig.  9).  Briefly 
summarizing  what  has  been  observed,  the  experimental  and  calculated  thickness 
profiles  agreed  very  well.  Where  there  was  a  deviation  between  theory  and 
experiment,  the  film  was  no  longer  SiC,  but  consisted  of  SiC  imbedded  in 
polycrystalline  silicon,  as  verified  by  X-ray  diffractometry.  Once  this  was 
accounted  for,  agreement  was  obtained  over  the  entire  deposition  range. 

A  number  of  SiC  depositions  were  made  and  the  highlights  are  listed  in 
Table  V.  The  structure  of  the  films  was  deduced  by  X-ray  diffractometry  and 
by  using  an  X-ray  Read  camera  which  is  similar  to  Debye-Scherrer/Laue  but 
optimized  for  thin  films.  The  X-ray  analysis,  resulting  in  what  is  called 
preferred  orientation,  produced  only  the  111  (and  higher  order)  reflections  of 
e-SiC.  The  refractive  index  of  the  films  was  measured  by  ellipsiometry.  Near 
the  stoichiometric  region  the  value  was  2.6,  agreeing  with  the  literature 
value. 

Film  composition  was  assessed  by  both  electron  microprobe  and  by  Auger 
analysis  and  good  agreement  was  found  for  the  two  techniques.  Composition 


profiles  for  Si  and  C  are  shown  in  Fig.  8.  By  varying  the  power  on  the  Si 
target,  the  stoichiometric  SiC  location  could  be  moved, as  indicated  in  the 
figure.  The  Auger  analysis  has  produced  noteworthy  results  on  the  Si  and  C 
line  shapes  and  how  the  line  shapes  are  influenced  by  the  various  environments 
of  the  Si  and  the  C  produced  by  the  variable-composition  depositions. 

For  the  SiC  deposited  on  AtgOj  .optical  measurements  could  be  made. 
Preliminary  measurements  of  the  transmitted  intensity  are  shown  in  Fig.  9. 

The  bandgap  evaluated  by  this  approximate  method  is  about  3  eV,  which  is  the 
bandgap  expected  for  SiC. 


Future  Direction 

Further  investigation  of  the  material  properties  of  the  described  SiC 
deposition  is  underway.  A  more  exacting  analysis  is  being  performed  to 
compute  the  bandgap  from  the  optical  data.  The  depositions  will  be 
scrutinized  by  SEM  and  TEM  to  establish  grain  size  and  growth  morphology. 
Electrical  propertis  of  the  film  will  be  measured,  particularly  resistivity 
and  carrier  type. 

As  mentioned  in  the  motivation  section,  a  new  set  of  experiments  is 
planned  to  investigate  a  surface  mobility-enhancement  effect.  First,  a  change 
from  dual  to  single  source  deposition  of  SiC  will  be  attempted  to  simplify  the 
evaluation  of  the  mobility-enhancement  effect.  The  single  source  is  powdered 
a-SiC  compressed  into  a  disk  shape.  Experimentally,  the  mobility  enhancement 
consists  of  co-sputtering  a  material  that  is  immiscible  with  both  the 
substrate  and  the  SiC  (Sn  is  a  candidate),  while  the  substrate  is  maintained 
above  the  melting  temperature  of  the  immiscible  material.  This  idea  is  shown 
schematically  for  a  dual  source  deposition  in  Fig.  10.  What  is  envisioned  is 
that  a  liquid-like  film  of  the  inert  material  will  form  a  2D  array  (liquid 
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lattice  idea)  just  above  the  substrate.  This  inert  array  will  enhance  the 
surface  mobility  of  the  sputtered  species  (Si  and  C)  (i)  by  perhaps  guiding 
the  sputtered  material  to  preferred  sites  for  epitaxial  film  growth,  or 
(ii)  by  being  in  the  liquid  state  and  thereby  allowing  higher  surface  mobility 
than  the  solid  substrate  would  allow.  The  effect  of  mobility  enhancement  has 
been  called  the  liquid-like  film  effect  (LLFE). 

A  mobility  enhancement  of  the  type  described  (the  mechanism  as  presented 
is  presently  fictitious)  has  been  seen  for  the  Cu-Pb  system  and  for  other 
systems  as  well.  Referring  to  Fig.  11,  as  the  Pb  vapor  pressure  is  increased, 
the  surface  diffusivity  (directly  related  to  the  surface  mobility)  of  copper 
is  increased.  Because  these  data  are  available,  another  set  of  experiments  is 
planned  in  parallel  with  the  SiC  +  Sn  depositions.  A  model  system  consisting 
of  a  simpler  deposition  reaction  will  be  investigated.  At  present,  two 
possible  model  systems  to  test  the  LLFE  are  Cu  sputtered  in  a  Pb  environment 
and  W  sputtered  in  a  Sn  environment. 

B.  Theoretical  Results 

In  present  atomistic  modeling  of  the  thermal  oxidation  of 
silicon,  (3-5),  structure  of  the  near  interface  region  is  of  crucial 
concern  because  this  is  a  key  determining  factor  in  the  formation  of  fixed 
oxide  charge  and  of  interface  charge  states  during  the  oxidation  process. 
Although  the  transition  of  Si  to  vitreous  SiOg  may  be  favored  to  proceed  via  a 
crystalline  a-cristobal ite  buffer  layer  containing  Si-interstitial 
defects^),  the  subsequent  internal  oxidation  of  these  Si j  defects  leads  to 
swelling  stresses  which  distort  the  lattice  structure  to  a  vitreous  structure 
of  a  very  similar  density.  An  understanding  of  this  important  transition  and 
of  the  residual  defects  stored  in  this  non-stoichiometric  buffer  layer 
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requires  an  understanding  of  the  differences  between  crystalline  and  vitreous 
Si02  on  an  atomistic  level.  Such  understanding  may  lead  to  ultimate 
controllable  reduction  in  the  magnitudes  of  the  fixed  oxide  charge,  Qf,  and  in 
the  density  of  Interface  charge  states,  N1t. 

The  thin  non-stoichiometrlc  interface  region  existing  during  thermal 
oxidation  is  expected  to  be  closely  analogous  to  bulk  Si 0X  films  formed  during 
reactive  sputtering  of  silicon  or  by  the  SIPOS  process.  Thus,  to  properly 
deal  with  this  class  of  structures,  one  needs  understanding  not  only  of  the 
crystalline  Sit^/vitreous  Si02  transition  but  also  of  the  structural 
consequences  associated  with  the  Introduction  of  SiOn  species  into  the  network 
and  Interstitial  spaces  of  the  S102 .  It  is  expected  that  this  will  aid  in 
understanding  the  roles  of  Cl,  H,  N  and  various  annealing  atmospheres  on  the 
structure  and  electrical  properties  of  SiOx  and  Si 02  films. 

To  understand  the  formation  and  properties  of  SiOx  films,  it  is  essential 
to  know  the  structure  of  these  films  at  an  atomic  level,  the  spatial  mechanism 
involved  in  the  formation  of  these  structures  and  the  energetics  of  these 
structures.  Since  direct  experimental  data  are  not  currently  available, 
recourse  must  be  made  to  theoretical  models.  Our  theoretical  approach 
Involves  two  steps:  (1)  modeling  of  vitreous  Si02  relative  to  crystalline 
Si02  and  (2)  modeling  vitreous  SiOx  structures  relative  to  vitreous  Si02.  The 
overall  problem  has  been  divided  into  a  sequence  of  simpler  tasks,  some  of 
which  have  been  completed  while  others  are  in  progress. 

(1)  Modeling  of  Amorphous  Silica 

The  current  models  of  vitreous  silica  are  not  sufficiently  detailed  for 
the  goals  envisioned,  thus,  a  different  approach  has  been  taken.  We  partition 
the  description  into  primary  building  block  considerations,  dimer  unit 
considerations  and  ring  structure  considerations.  Via  this  path,  one  is  able 
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to  generate  any  possible  macroscopic  structure  be  it  crystalline  or  vitreous, 
be  it  stoichiometric  or  non-stolchiometric  and  one  is  more  able  to  understand 
the  atomic  scale  defects. 

(a)  Basic  Units 

All  silica  structures,  crystalline  or  vitreous,  are  made  up  of  the 
same  basic  units  -  oxygen  tetrahedra  with  a  silicon  Interstitial.  Adjacent 
tetrahedra  are  linked  by  the  sharing  of  a  corner  referred  to  as  the  bridging 
oxygen.  In  perfect  silica,  all  the  corners  of  the  tetrahedra  are  bridging 
oxygens.  In  real  silica,  some  of  these  bridging  oxygens  will  have  moved  into 
interstitial  positions  leaving  vacancies  at  these  sites.  Thus,  the 
description  of  any  silica  structure  essentially  involves  describing  the 
relative  orientation  of  every  pair  of  adjacent  linked  tetrahedra  which  we  will 
refer  to  as  the  dimer  unit. 

The  tetrahedra  are  the  least  variable  units  In  the  silica  structures  and 
practically  retain  their  shape  and  dimensions  In  most  silicas,  even  when 
subject  to  large  pressures.  This  view  Is  reinforced  by  noting  that  the  free 
energy  of  formation  of  SiC^  at  room  temperature  is  212  kcal/mole  whereas  the 
free  energy  difference  between  any  two  silica  structures  is  only 
~  2-4  kcal/mole^6).  Thus,  ~99%  of  the  free  energy  storage  resides  in  the 
tetrahedral  units. 

The  tetrahedra  are  almost  perfectly  regular  as  Illustrated  in  Fig.  12. 
The  six  edges  are  nearly  equal  (2.65  ±  0.10  A);  the  four  corners  are  at 
nearly  the  same  distance  from  the  central  Si  (1.62  t  0.05  A).  Thus,  the  0- 
S1-0  angles  are  about  110  t  10°.  The  small  deviations  from  regularity  and 
the  small  variations  in  shape  and  dimensions  of  the  tetrahedra  in  different 
silica  structures  should  be  considered  as  second  order  effects.  These  are 
evidently  due  to  the  different  environments  of  the  unit  in  different 
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directions  and  to  the  non-uniform  forces  present.  We  will  ignore  these  small 
variations  in  our  initial  considerations. 

A  scheme  for  linking  the  tetrahedra  into  a  macroscopic  structure  is 
illustrated  in  Fig.  13.  We  start  with  a  single  tetrahedron  (#1)  and  we  link 
it  to  four  others  as  first  nearest  neighbours  (2,3,4  and  5)  by  defining  the 
three  angular  parameters  specifying  each  dimer  (1-2;  1-3;  etc.).  Each  of 
these  four  unis  can  then  be  linked  to  three  new  mits  leading  to  12  second 
nearest  neighbours.  This  procedure  can  be  continued  until  an  adequate  volume 
of  the  structure  has  been  built.  In  order  for  this  procedure  to  be  manageable 
and  meaningful,  restrictive  criteria  and  classifying  procedures  are 
necessary.  There  is  no  general  reason  to  expect  the  environment  of  any  unit 
in  any  given  structure  to  be  vastly  different  from  that  of  another  in  the  same 
structure.  Thus,  it  should  be  adequate  to  define  only  a  few  dimers  in  a 
sample  volume.  For  most  crystalline  silica  structures,  defining  a  single 
dimer  unit  and  a  few  symmetry  rules  has  been  shown  to  be  sufficient  to  build 
the  entire  structure  to  a  good  approximation 

(b)  Dimer  Unit 

Figure  14  illustrates  how  the  three  angles  (e,«lt62)  defining  a 
dimer  unit  can  be  chosen.  The  angle  8  is  made  by  the  two  silicons  at  the 
bridging  oxygen.  The  angles  5^,  «2  are  the  tw0  dihedral  angles  defining  the 
position  of  one  oxygen  from  each  of  the  two  tetrahedra,  0^  and  02 
respectively,  relative  to  the  Si j-Oq-SI 2  plane.  Since  the  size  and  shape  of 
the  tetrahedra  are  known,  the  other  four  oxygens  are  automatically  fixed  and 
the  dimer  Is  fully  defined  by  4,  6 ^  and  6^  .  The  dihedral  angles  defining 
O3  and  O4  can  be  expected  to  be  about  5^  ±  120°.  Similarly  05  and 
0g  are  defined  by  dihedral  angles  ■  fi2  4  120°. 
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The  bridge  angle  8  is  generally  found  to  be  close  to  144u  in  most 
crystalline  silica^8).  The  current  models  of  amorphous  Si02  are  based  on  wire 
models,  sometimes  relaxed  by  computer  simulation,  which  assume  a 
distribution  of  0  from  120°  to  180u  with  a  maximum  at  144°  and  random 
values  of  Sij  and  Si 2-  However,  the  range  of  0  in  vitreous  Si02  is  probably 
much  less  than  the  assumed  range  of  120u-180°.  In  very  early  data  on  high 
cristobalite,  0  was  assumed  to  be  180u  but  recent  determinations  of 
e  -  144°  Many  investigators  still  use  the  old  incorrect  value.  Data 

on  high  tridymite  show  a  similar  shift  in  8  on  redetermination.  It  has 
recently  been  suggested  that  non- bonded  repulsion  between  the  silicons  in 
a  dimer  determines  their  distance  of  closest  approach.  For  a  Si-0  distance  of 
~  1.62  A,  the  Si-Si  separation  is  ~  3.08  A  for  0  =  144°.  The  Si-Si 
repulsion  becomes  severe  when  0  is  reduced  to  ~  120u  (Si-Si  separation 
~  2.8  A)  and  0-values  lower  than  this  are  thus  practically  ruled  out. 

The  restrictions  on  the  6's  will  be  essentially  due  to  interactions 
between  the  neighboring  oxygens.  However,  the  closest  separation  observed 
between  oxygens  in  most  Si02  structures  is  much  larger  than  the  nonbonded 
contact  separation  between  oxygens  in  a  tetrahedron.  Thus,  restrictions  due 
to  non-bonded  repulsions  between  oxygens  is  not  likely  to  be  an  important 
factor  in  silica  structure. 


(c)  Ring  Structures 

At  this  level  of  structural  organization,  we  begin  to  see  a 
possible  multiplicity  of  rings  constructed  from  the  basic  and  dimer  units.  To 
build  an  N-membered  oxygen  ring,  we  start  with  tetrahedral  unit  1,  link  it  to 
unit  2  and  link  unit  2  to  unit  3,  etc.,  as  in  Fig.  15.  If  the  (N  +  l)th  unit 
is  the  same  as  the  first,  then  we  have  a  closed  ring  of  N  bridging  oxygens. 


If  it  is  not  the  same,  then  we  have  an  open  ring  of  N  bridging  oxygens.  These 
rings  have  approximately  equal  sides  which  are  equal  to  the  edges  of  the  basic 
tetrahedral  unit  (~  2.65  A)  and  are  a  characteristic  feature  of  silica 
structures,  both  crystalline  and  amorphous.  In  a  sense,  these  rings  sample 
the  structure  over  a  small  domain  and  it  is  useful  to  describe  silica 
structures  in  terms  of  the  oxygen  rings. 

It  is  obvious  that  there  should  be  consideable  restrictions  on  the  6's  in 
order  to  ensure  ring  »:lojjre  and  to  ensure  that  the  e  for  the  dimers  of  the 
ring  are  in  the  acceptable  range.  We  have  developed  effective  procedures  by 
which  one  can  identify  ’.he  ranges  in  6  that  lead  to  the  formation  of  closed 
rings  of  differ©';.  *es  and  different  geometries;  i.e.,  containing  3,4,5  or  6 
bridging  oxygens.  Although  the  three  angular  parameters  9,  and  were 
sufficient  to  describe  the  dimers,  for  the  rings,  one  needs  two  sets  of 
parameters:  (1)  those  defining  the  geometry  of  the  oxygen  ring  and  (2)  those 
defining  the  orientation  of  the  constituent  tetrahedra  relative  to  the  ring. 

We  will  first  illustrate  the  use  of  these  parameters  with  a  6-ring. 

(1)  The  oxygens  of  a  6-ring  are  identified  as  1,2 . 6  and  the  sides 

1-2,  2-3 . 5-6,  6-1  are  all  roughly  equal  since  they  are  an  edge  of  a 

tetrahedron,  the  six  angles  at  the  vertices  1,2, . 6  of  the  ring  or, 

equivalently,  the  separation  distances  between  alternate  oxygens  are  needed 
parameters.  These  separations  are  between  oxygens  from  adjacent  tetrahedra 
forming  a  dimer.  Three  dihedral  angles  complete  the  set  of  parameters 
required  to  define  the  ring  geometry,  these  dihedral  angles  could  be 
equivalently  described  by  the  separations  1-4;  2-5  and  3-6.  The  geometry  of 
the  ring  in  o-quartz  at  a  pressure  of  68  Pa  is  shown  in  Fig.  16. 

(2)  The  geometry  of  the  ring,  as  define  above,  fixes  only  one  edge  of 
the  six  constituent  tetrahedra.  Assuming  the  tetrahedra  to  be  regular,  we 


need  only  one  other  parameter  to  define  the  tetrahedra;  i.e.,  a  dihedral  angle 
around  the  concerned  edge  with  a  chosen  reference  plane  in  the  ring.  This 
dihedral  angle,  which  we  shall  call  S',  can  be  used  to  define  the 
corresponding  Si  which  will  automatically  lead  to  a  definition  of  the  two 
remaining  oxygens  in  the  tetrahedron. 

Using  a  different  numbering  scheme  to  identify  the  oxygens  in  Fig.  17,  we 
will  illustrate  the  procedure  being  used  for  classifying  the  ring  structure  of 
size  3,4,5  and  6.  By  varying  the  angles  e,6j,62  =  $,6p62  for  the  dimer 
(Fig.  17a)  over  say  0  -  120°  -  180°,  6j  -  0-180°  and  «2  -  0-360°,  all  the  0-0 
distances  may  be  calculated  between  the  two  tetrahedra  forming  the  dimer.  In 
the  complementary  notation,  the  angle  ^  at  1  is  varied  from  30°  to  160°. 

By  using  dihedral  angles  representing  rotation  about  1-2  and  1-3,  with  respect 
to  plane  312  as  reference,  the  silicons  are  located  and  the  Si-Si  distance  in 
the  dimer  found.  The  dimer  is  now  completey  defined. 

For  each  ^  ,  the  accepable  Si-Si  separation  will  place  varying  amounts 
of  restriction  on  the  6 ' ' s  and  the  distance  2-3  is  known.  Clearly,  when  2-3 
equals  1-2,  we  have  a  3-oxygen  ring.  Considering  Fig.  17b,  the  oxygen  4  is 
defined  by  $2  and  a  dihedral  angle  (3124).  Similarly,  oxygen  5  is  defined  by 
and  the  dihedral  angle  (2135).  When  the  distance  3-4  or  2-5  is  zero,  we 
have  a  3-oxygen  ring.  When  the  distance  of  either  of  them  is  equal  to  the  1-2 
distance,  we  have  a  4-oxygen  ring.  When  the  distance  4-5  =  distance  1-2,  we 
have  a  5-oxygen  ring.  When  4-5  =  2-3,  we  have  a  6-oxygen  ring  and  0g  is 
defined  by  another  dihedral  angle  around  5-4.  In  each  case,  our  earlier  study 
of  the  dimer  will  restrict  the  dihedral  angle  values  because  of  the 
restriction  of  the  range  of  the  bridge  angle  (or  of  the  Si-Si  separation). 

Although  the  results  of  these  studies  are  not  yet  fully  analyzed,  we  will 
briefly  describe  the  trends  that  are  appearing:  (i)  a  3-ring  is  highly 
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restricted  in  regard  to  all  of  its  parameters.  The  ring  angles  should  all  be 
close  to  60u  since  all  the  sides  are  nearly  equal,  the  silicons  should  be 
located  almost  on  the  plane  of  the  ring  and  outside  of  the  ring  to  ensure 
adequate  Si-Si  separation  (see  Fig.  18).  The  bridge  angles  are  about  130°  in 
this  case  which  is  probably  a  strained  configuration. 

(ii)  A  4-ring  has  a  little  more  flexibility  both  with  respect  to 

and  the  6 '  ‘ s .  The  restriction  on  4  and  6'  ,  are  even  more  relaxed  in  the 
cases  of  5-rings  and  6-rings.  All  such  restrictions  will  be  even  more  severe 
if  the  acceptable  range  in  the  bridge  angles  is  less  than  the  assumed  range 
(120°-180°]. 

(iii)  It  appears  that  dimer  parameters  leading  to  rings  other  than  6  are 
unlikely  to  exhibit  long  range  order.  However,  formation  of  6-rings  does  not 
necessarily  imply  long  range  order  although  most  crystalline  Si O2  can  be 
described  in  terms  of  6-rings^1)  with  the  ring  geometry  differing  for 
different  polymorphs. 

(iv)  The  density  of  the  structure  is  another  factor  in  its 
acceptability.  Rings  like  those  in  o-cristobal ite  may  be  found  in  the 
vitreous  Si02  structure  because  the  densities  are  similar  whereas,  rings  close 
to  those  found  in  quartz  can  be  ruled  out  for  vitreous  Si02  because  of  the 
much  higher  density  for  quartz. 

(v)  One  other  criterion  still  to  be  considered  for  the  acceptibil ity  of 
a  stucture  is  the  0-0  separation.  As  mentioned  earlier,  all  oxygens  in 
crystalline  silica  have  large  enough  separations  that  non-bonded  repulsions 
are  not  a  factor.  The  repulsive  forces  due  to  residual  charges  on  the  0's 
will  ultimately  decide  the  resultant  geometry  of  the  real  structures.  An 
analysis  of  the  neighboring  oxygen  distances  in  crystalline  Si02  will  provide 
some  clues  in  this  regard.  For  example,  the  structure  of  a-quartz  has  been 


determined  at  different  pressures  from  0.001  Pa  to  68.0  Pa^).  The  results 
indicate  that  the  tetrahedra  suffer  very  little  deformation  and  the 
compression  is  brought  about  mostly  by  the  bending  of  the  bridge  angle, 9. 


Over  this  pressure  range,  the  Si-0  distance  changes  from  1.615  A  to 
1.601  A;  the  Si-Si  separation  drops  from  3.06  A  to  2.95  A  and  e  changes  from 
142°  to  135°  .  Looking  at  the  smallest  0-0  distance  in  Fig.  5,  it  drops  from 
3.30  A  to  2.84  A  and  it  is  seen  that  the  closest  oxygens  are  from  second 
neighbor  tetrahedra. 

(2)  Modeling  the  Energetic  Aspect  of  Interstitial  Species  in  silica 

In  the  present  modeling  of  thermal  oxidation  of  silicon,  the  emphasis 
has  shifted  from  a  parametric  description  to  an  atomistic  description.  At 
this  new  level  of  modeling  we  are  confronted  with  a  need  for  the  following 
information:  (i)  solubility  data  for  the  various  species,  (ii)  dissociation 
and  ionization  data  for  the  various  species  within  the  Si02,  (iii)  bridging 
oxygen  vacancy,  VQ,  and  Frenkel  defect,  V0  +  0^,  formation  energies, 

(iv)  reaction  energies  for  the  internal  oxidation  of  Si ^  species  to  SiO, 

Si02,  Si03  and  (v)  binding  energies  for  the  V0  -  Xx  pairs.  In  addition  to 
needing  this  information  for  the  bulk  Si02,  it  is  also  needed  for  the  non- 
stoichiometric  interface  blocking  layer  of  Si02  which  is  thought  to  be  so 
critical  to  the  velocity  of  oxidation  and  to  the  formation  of  the  fixed  oxide 
charge,  Qf.  This  foregoing  information  is  essential  to  an  understanding  of 
the  mechanisms  for  (a)  formation,  (b)  Ct  adsorption,  (c)  Qf  annealing,  and 
(d)  network  versus  interstitial  diffusion  for  various  species.  In  addition, 
it  is  needed  in  order  to  determine  the  population  of  ionized  versus  neutral 
species  so  that  a  meaningful  assessment  may  be  made  of  electrical  conduction 
in  Si02  at  oxidizing  temperatures  as  well  as  of  electric  field  effects  on 


oxidation.  Some  portion  of  this  information  is  accessable  by  experimentation; 
however,  most  of  it  is  not  and  the  only  recourse  we  have  is  to  theoretical 
computation.  It  is  for  this  reason  that  we  have  used  a  molecular  orbital 
approach  to  resolve  these  questions. 

We  have  used  the  CNDO/2  method  because  the  computation  costs  are 
considerably  smaller  than  ab  initio  results^1  .  Some  limitations  of  our 
present  method  are  that  (i)  the  Si-0  bond  lengths  are  not  optimized,  for 
minimum  total  energies,  (ii)  the  reaction  energies  are  determined  at  0°K  , 

(iii)  these  calculated  energies  are  not  absolute  so  that  only  relative 
energies  between  different  species  or  conditions  are  likely  to  be  relevant, 

(iv)  we  can  only  handle  82  orbitals  and  (v)  we  can  only  handle  closed  shell 
configurations.  In  our  present  method,  we  have  adopted  de  Jong's  ^  criteria 
that  (1)  all  orbital  energies  of  occupied  MO's  should  be  negative,  (2)  the 
total  energy,  Ey,  should  be  minimized  at  a  value  of  e  -  140°  ,  (3)  the 
binding  energy  of  the  molecule  should  be  negative  or  close  to  zero  and 

(4)  bond  overlap  populations  and  other  measures  of  bond  energies  should  be 
consistent  with  geometries  and  dimensions  observed  in  silicate  structures. 

In  the  previous  report  we  achieved  the  following  results  for  0  and 
Si  interaction  with  a  dimer  unit  and  a  tetramer  ring: 


0  +  H6Si207  Z  HgSi 20g  ;  4Ey  =  -6.625  eV  (1) 

0  +  2e"  +  HgSi207  J  HgSigOg2  ;  AEy  =  +  6.76  eV  (2) 

Si  +  HgSi207  j  HgSi307  ;  4Ey  =  +  8.717  eV  (3) 

o+ 

Si  +  HgSi207  ;  HgSi307  +  2e  ;  4Ey  =  +  19.2  eV  (4) 
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In  the  above  equations,  the  electron  work  function  was  chosen  as  5.1  eV.  A 
larger  value,  by  6E,  would  increase  AEy  in  Eqs.  (2)  and  (4)  by  25E.  The 
following  results  with  the  tetramer  ring  (HgSi^Q^)  were  also  achieved  (with 
0.j  at  x  =  0.5  A) 

0  +  tetramer  +■  tetramer  +  0(x  =  0.5);  aE-j.  =  -3.042  eV  (5) 

0  +  2e~  +  tetramer  t  tetramer  +  Cf  (x  =  0.5);  AEy  =  +9.9  eV  (6) 

0  +  tetramer  X  tetramer  0++  +  2e"(x  =  0.5);  AEy  =  +9.1  eV  (7) 

Again,  the  *<ork  function  was  chosen  as  4.5  eV  in  Eqs.  (6)  and  (7). 

From  Eqs.  (1-7)  we  note  that,  if  monatomic  oxygen  is  available  in  the  gas 

phase,  it  will  be  strongly  absorbed  by  both  the  dimer  and  the  tetramer.  This 

conclusion  also  applies  to  0=  (add  -  9  eV  to  Eqs.  (2)  and  (6)).  However, 

ionization  insitu  to  form  0=  is  not  a  high  probability  event.  The  difference 

in  energetics  between  the  dimer  calculation  and  the  tetramer  calculation  is 

probably  the  same  as  that  found  by  placing  the  0  above  versus  below  the 

bridging  oxygen  in  the  dimer.  The  reduced  approach  angle  provides  the 

enhanced  repulsion  force.  This  proposal  will  be  tested  shortly.  Thus,  the 

AEy  values  given  by  Eqs.  (1)  and  (5)  apply  to  placing  the  0  adjacent  to  the 

obtuse  angle  side  versus  the  acute  angle  side  respectively  of  a  bridging 

oxygen  in  a  ring  interior  to  the  Si02  volume.  It  is  also  very  likely  that  the 

difference  in  these  two  values  of  AEy  is  close  to  the  activation  energy  for 

diffusion  of  the  0  species  through  a  Si structure  made  from  4-membered 
o 

rings;  i.e.,  Q  -  3.5  eV,  via  the  interstitial  spaces. 
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The  main  question  addressed  during  the  present  reporting  period  involved 
the  interaction  of  O2  and  i^O  with  the  Si02  structure  and  whether  or  not 
dissociation  is  probable.  First,  we  considered  the  C^-tetramer  reaction  and 
found  that  the  minimum  energy  condition  involved  the  O2  being  placed  at  the 
center  of  the  ring  with  the  0-0  bond  being  perpendicular  to  the  plane  of  the 
ring  and  fixing  its  length  to  be  that  in  O2  at  25°C .  The  following  energetic 
reaction  was  observed 

0 2  +  tetramer  Z  (tetramer  +  0  2in  center);  AE^  =  +  7.387  eV  (8) 

This,  obviously,  is  not  the  energetically  most  favorable  site  in  the  Si O2 
volume,  even  though  it  is  in  the  4-ring,  because  it  would  lead  to  an 
impossibly  small  solid  solubility  for  the  02  compared  to  the  experimental 
value  of  ~  1016per  cc  at  1000°C  (15).  We  shall  see  later  that  a  much  more 

favored  site  occurs  on  the  outside  of  the  4-ring  where  the  02  can  approach  the 

obtuse  angle  side  of  the  bridging  oxygen.  For  the  moment,  let  us  consider  the 

dissociation  of  the  O2  in  Eq.  (8)  to  form  20,  one  of  which  enters  an  adjacent 

tetramer  and  both  0's  are  located  at  0.5  A  from  the  cente*-  a?  :ts  tetrad?-; 

i  .e. , 


02  +  2  tetramers  t  2(tetramer  +  0(x=0.5));  AEj  =  +2.827  eV  (9) 

With  this  reaction  we  can  expect  a  much  greater  solubility  for  02  in  Si O2 
although  AEy  is  still  too  large.  Perhaps  with  a  6-ring,  the  value  of 
AET  would  be  sufficiently  reduced  to  match  the  experimental  data  (when  entropy 
factors  are  also  included). 

When  we  turn  to  the  interaction  of  O2  with  the  dimers,  a  very  interesting 
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result  is  obtained  depending  upon  the  orientation  of  the  02  relative  to  the 
plane  of  the  Si -0-Si  atoms  in  the  dimer.  These  results  are  illustrated  in 
Fig.  19  where  Ey  is  plotted  as  a  function  of  the  distance  z  between  the 
bridging  oxygen  and  the  0-0  line  for  two  orientations  of  the  02  molecule  (a) 
in  the  Si -0-Si  plane  of  the  dimer  and  (b)  perpendicular  to  the  Si -0-Si  plane 
of  the  dimer.  In  the  former  case,  we  have  the  reaction 

02(p)  +  HgSi 207  j  HgSi209  ;  AEy  -+  0.80  eV  (10) 

and,  from  Fig.  1,  we  see  that  as  Z  is  reduced  from  large  values,  there  is  a 
slight  dip  in  AEy  and  then  an  increase  as  z  is  reduced  to  ~  1.21  A  which  is 
close  to  dQ_g  =  1.2074  A  .  Turning  to  case  (b),  we  note  that  if  the  02(p) 

molecule  rotates  90u  to  the  02(b)  position  at  z  ~  1.95  A,  a  second  branch  is 

available  for  Ey  and  this  is  a  steeply  descending  branch  (see  Fig.  19).  As  z 
decreases  and  Ey  descends  along  branch  (b),  the  magnitude  of  the  Si-0  bridge 
bond  decreases  and  reaches  a  zero  value  at  z  <  1.36  A  and  a  minimum  at  z  = 
1.15  A  .  For  this  case  the  reaction  equation  is 

02(jl)  +  H6Si207  t  H6si2°9;  aET  “  '  3,54  eV  (11) 

This  is  a  bond  breaking  reaction  which  on  the  surface,  appears  to  be  a 
mechanism  for  vacancy  formation  at  the  bridging  oxygen  site.  However,  before 
we  can  be  confident  of  the  reality  of  this  result,  we  must  ascertain  whether 
it  is  possible  for  the  02  molecule  to  rotate  by  90u  to  form  the  02(i) 
configuration  within  the  volume  of  the  Si02  or  on  the  surface  of  the  Si 02 

because  of  the  hindrance  of  the  basic  tetrahedral  units  of  the  dimer  at  their 

normal  dihedral  angles,  6^  and  5^.  Although  it  may  be  unlikely  for  this 


normal  dihedral  angles,  and  6g.  Although  it  may  be  unlikely  for  this 

interesting  event  to  take  place  in  the  volume  of  Si02,  it  is  likely  to  occur 
on  the  surface  of  Si02 

In  Fig.  20,  the  Ey  values  are  given  for  the  interaction  of  an  H20 
molecule  with  the  dimer  unit  for  two  configurations  of  the  H20  using 
<  HOH  =150°  and  dg_H  *  0.96  A.  Only  the  (b)  branch  is  a  bond  forming 
reaction.  The  following  reaction  equations  apply  to  these  two  branches 

HgO(a)  +  HgSi207  £  H8Si20g(a);  aEt  =  +  1.42  eV  (12) 

HgO(b)  +  HgSi g07  *  HgSi20g(b);  AEy  =  -0.17  eV  (13) 

A  third  configuration  (c)  where  a  single  H  was  bonded  to  the  bridging  oxygen 
was  initially  investigated  by  De  Jong^7)  who  found 

H20(c)  +  H6Si207  1  HgSi20g(c)  ;  AEy  =  -0.35  eV  (14) 

Once  again,  we  would  expect  these  three  configurations  to  fit  on  the  Si 02 
surface  but  there  are  real  doubts  that  they  would  fit  in  the  bulk  Si02 
without  some  repulsive  energy  contributions. 

An  additional  reaction  of  H20  with  Si02  that  can  be  simply  assessed  with 
the  dimer  is  the  dissociation  reaction  into  H+  and  OH";  i.e., 

H20  +  2HgSi207  +  HgSigOy (H+)  +  HgSigO^OH") ;  AEy  =  +5.18  eV  (15) 

In  Eq.  (15),  the  H+  and  the  OH"  are  absorbed  at  the  bridging  oxygen  sites  at 
their  minimum  energy.  This  large  positive  energy  change  arises  from  the  large 


dissociation  energy  of  H20  into  H+  and  OH"  since  de  Jong^7)  found  the 
following  reactions 

OH”  +  HgSi207  t  H6Si207(OH”);  AEy  =  -12.65  eV  (16) 

H+  +  HgSi207  t  H6Si207(H+);  AEy  =  -0.65  eV  (17) 

It  is  the  large  difference  in  the  binding  energies  of  OH-  vs  H+  that  leads  to 
the  strong  negative  surface  charge  of  silica  associated  with  the  isoelectric 
point  of  ~  2.0. 

If  we  next  consider  the  possible  reactions  of  H20  (-19.87  a.u.)  with  the 
tetramer  unit  (-236.5368  a.u.),  we  find  the  most  favorable  reaction  is  that 
given  by  configuration  (g)  in  Fig.  21.  Here,  the  0  is  in  the  center  of  the 
ring  while  the  H's  are  above  and  below  the  ring  directed  towards  an  0  unit  of 
the  tetramer  ring.  The  reaction  equation  is 

H20(g)  +  tetramer  £  tetramer  (+H20(g));  AEy  =  +9.84  eV  (18) 

We  also  have  the  dissociation  and  ion  binding  reactions  as 
H20  +  2  tetramers  J  tetramer(H+)  +  tetramer(OH’) ;  AEy  =  +16.28  eV  (19) 


H  +  tetramer  £  tetramer  (H+)  +  e”;  AEy  =  -10.28  eV  (20) 

OH’  +  tetramer  t  tetramer  (0H‘)  ;  AEy  =  +  1.88  eV  (21) 
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Once  again,  the  large  difference  between  Eq.  (18)  and  the  sum  of  Eqs.  (19)  and 
(20)  is  the  dissociation  energy  of  H20  into  H  and  OH.  In  Eq.  (20),  H  is  taken 
to  have  zero  energy  by  convention. 

From  the  foregoing,  we  can  note  that  the  calculations  involving 
interactions  with  the  dimer  unit  are  (a)  completely  appropriate  for  Si02 
surface  in  contact  with  the  gas,  (b)  approximately  correct  for  the  bulk  when 
we  have  a  small  interstitial  species  interacting  with  the  obtuse  angle  of  the 
bridging  oxygen  and  (c)  incorrect  for  the  Si 02  bulk  when  we  have  a  large 
interstitial  species  interacting  with  the  obtuse  angle  side  of  the  bridging 
oxygen  because  of  neglected  interactions  with  the  surrounding  environment. 

The  interactions  with  the  tetramer  ring  give  some  insight  into  normal  bulk 
Si02  interactions;  however,  we  really  need  to  treat  the  actual  6-ring  to 
obtain  final  numbers.  The  tetramer  calculations  serve  as  a  valuable  and 
inexpensive  learning  tool.  Tabulated  formation  energies  for  species 
interaction  with  the  dimer  and  tetramer  units  are  given  in  Table  VI. 

In  closing  this  section,  a  word  may  be  said  concerning  Si 02  spectral 
information  generated  by  the  foregoing  type  of  calculations.  Although  CNDO/2 
calculations  do  not  give  the  0(1S)  core  level  energy,  shifts  in  this  energy 
level  may  be  inferred  from  changes  in  the  calculated  gross  atomic  population 
of  oxygen ( 16 K  For  example,  from  the  calculated  nearly  neutral  charge  on  the 
02  in  the  center  of  the  tetramer,  it  can  be  inferred  that  the  0 ( IS)  spectrum 
from  interstitial  02  should  be  almost  identical  with  that  of  molecular  02. 
However,  from  the  calculated  charge  on  the  bridging  oxygen,  the  0 ( IS)  core 
level  shifts  should  be  less  than  ~  0.5  eV  .  Thus,  a  large  solubility  of  02 
molecules  in  Si02  may  produce  a  shift  to  more  negative  energies  (~  13  eV)  for 
the  0 ( IS)  photoemission  spectrum. 


-27- 


TABLE  I 


SUMMARY  OF  EXPERIMENTS 
ON 

REACTIVE  DEPOSITION  OF  ZnOx,  ZrOx>  AiOx 


Oxide  Systems 

Structure 

AHd(kcal 

Si/O  yes 

Amorphous 

210 

(ZnO)  Zn/O  yes 

Xtal  (Hex) 

166.8 

(Zr02)  Zr/O  yes 

Xtal  (Cubic) 

259.5 

a  -  AX./0  yes 

Amorphous 

266.8 

*  -  VA("°x> 

<P  " 

V"> 

AHd  -  Heat  of  decomposition  (kcal) 
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TABLE  II 


SUMMARY  OF  EXPERIMENTS  ON  ASl Ox  SYNTHESIS 


A* 

Ar(torr) 

02(torr) 

Comments 

194 

500  co 

(5.8  A/sec) 

4  x  i(f3 

2  x  io'4 

At,  At0x 

195 

250 

(3.2  A/sec) 

8  x  io'3 

1.6  x  10'4 

Oxidized  target 

196 

400 

(4-1  A/sec) 

8  x  io'3 

2  x  10‘4 

Oxide  (clear) 

197 

650 

(7.0  A/sec) 

8  x  10'3 

2.2  x  io‘4 

Oxide,  AtOx 

198 

650 

(7.0  A/sec) 

8  x  io"3 

2.5  x  10"4 

-A 

Oxide,  AtOx 

At  -  41.7 

199 

650 

(7.0  A/sec) 

8  x  10  J 

0/A  =  1.40, 

3.5  x  10  4 
nf  =  20 

AtO 

x  0  -  58.3 
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TABLE  III 


EXPERIMENTAL  PARAMETERS  USED  IN  THE  SYNTHESIS  OF  ZrOx 


Run 

Zr( Power) 

02(torr) 

Ar(torr) 

Comments 

029 

700  watts 

** 

1 

O 

r— 

X 

CVI 

8  x  10‘3 

Zr°2"Cublc 

030 

700 

i 

o 

X 

8  x  lo"3 

Zr°2"Cublc 

031 

700 

7  x  10‘5 

8  x  IQ"3 

Zr°2"Cubic 
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TABLE  IV 


EXPERIMENTAL 

PARAMETERS  USED  IN 

THE  SYNTHESIS 

OF  ZnOx 

Run 

Z  (Pwr) 
n 

OgUorr) 

Ar(torr) 

Conments 

022 

150 

6.5  x  10'4 

8  x 

I0'3 

target  oxidized 

023 

300 

6.5  x  10  4 

8  x 

10‘3 

reduced  ZnOx 

024 

250 

5.7  x  10"4 

8  x 

10"3 

reduced  ZnOx 

025 

250 

6.5  x  10'4 

8  x 

10‘3 

ZnO 

026 

250 

6.5  x  10"4 

8  x 

10'3 

ZnO 
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TABLE  V 

HIGHLIGHTS  OF  SiC  DEPOSITIONS 


750°C 

830*  C 


SiC/Si 

SiC/A£203 

SiC/Si 

SiC/Si 


amorphous 


B  SiC  111  reflection 
"preferred  orientation" 

B  SiC  111  reflection 


B  SiC  111  reflection 
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TABLE  VI 
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INTERACTION  ENERGIES  FOR  INTERSTITIAL  SPECIES  INTERACTIONS 
WITH  DIMER  (D)  AND  TETRAMER  (T) 


Species 

Dimer  AET° ( eV ) 

Tetramer  AE-^  (eV) 

Si0 

+  8.72 

Si^ 

+19.2 

0 

-  6.63 

-  3.04 

0= 

+  6.76 

+  9.9 

_*f+ 

0 

+  9.1 

02(p) 

+  0.80 

0(x) 

-  3.54 

+  7.4 

02  +  2T  t  2(T+0) 

+  2.83 

H20(a) 

+  1.42 

H20(b) 

-  0.17 

H20(c) 

-  0.35 

H20(g) 

+  9.84 

H+ 

-  0.65 

-10.28 

OH' 

-12.65 

+  1.88 

H20  +  2T  *  TH+  +  TOH* 

+  5.18 

+16.28 
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FIGURE  CAPTIONS 


1.  (a)  Schematic  diagram  of  new  sputtering  chamber,  (b)  Schematic  of 
new  gas  dispersion  ring. 

2.  Comparison  of  our  refractive  index  vs.  stoichiometry  data  with  the 
SIPOS  data. 

3.  Controlling  Pn  values  to  produce  a  given  film  stoichiometry  (x  in 

u2 

SiOx)  for  several  sputtering  power  levels. 

4.  Variation  of  index  of  refraction,  n»  and  'film  thickness,  t^,  as  a 
function  of  radial  distance  from  the  edge  of  the  wafer  for  a  Si0-|.5 
film. 

5.  Plot  of  0/A£  ratio  in  A£0x  film  as  a  function  of  radial  distance  from 

the  center  of  the  wafer  for  sample  197  (A£  -  650  watts,  02  -  2.2  x  10’4 
-3  0 

torr,  Ar  -  8  x  10  torr,  and  7.0  A/sec  deposition  rate). 

6.  Plot  of  electrical  resistance  as  a  function  of  radial  distance  for 
ZnOx  film  (sample  81-024). 

7.  Schematic  illustration  of  our  two-source  deposition  system  for  Si Cx 
film  formation  (SC  silicon  and  graphite  sources,  4  micron  Ar  pressure, 
25°C  -  830°C  substrate  temperature,  and  (100)  Si  or  (1102)  A£203 
substrates) . 

8.  Variation  of  C  and  Si  concentration,  evaluated  via  electron  microscopy, 
with  distance  along  the  SiCx  film  for  two  different  deposition  power 
levels  for  the  Si  (250  W  and  500  W). 

9.  Log  transmitted  intensity  as  a  function  of  wavelength  to  measure  the 
SiC  bandgap. 

10.  Schematic  illustration  of  Sn  liquid-like  film  effect  on  the  deposition 
of  SiC. 
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11.  Temperature  dependence  of  surface  diffusion  coefficient  of  Cu,  D$,  at 
several  different  values  of  Pb  vapor  pressure.  (Adapted  from  Fig.  16 
and  Fig.  26  of  Ref.  2. ) 

12.  Basic  structural  unit  in  Si02:  oxygen  tretrahedron  plus  silicon 

interstitial.  The  tetrahedron  is  almost  regular  (4  Si-0  distances  = 

0  0 
1.62  ±  0.05  A;  6  0-0  distances  a  2.65  ±  0.10  A,  and  0-Si -0  angles  « 

110°  ±  10°). 

13.  Scheme  for  linking  tetrahedra:  1  is  the  starting  unit,  2-5  are  first 
neighbors  of  1  and  6-17  are  second  neighbors  of  1.  Three  angular 
parameters  define  the  relative  orientations  of  the  linking  tetrahedra. 

14.  (a)  Drawing  of  a  tetrahedra  dimer:  angles  e,  6-j  and  62  fully  define 
the  dimer,  (b)  The  dihedral  angle,  61 ,  (the  distance  Si2  -  0-j  is  a 
measure  of 

15.  Scheme  used  for  the  formation  of  closed  rings. 

16.  Illustration  of  key  angles  and  distances  in  the  6-ring  of  a-quartz 

0  0 

at  68  P^  (distance  2-5  is  4.89  A  and  3-6  is  4.94  A).  At  0.001  P^, 

0 

the  distance  1-4  is  3.30  A. 

i  i 

17.  (a)  Classification  of  the  dimer  structure  via  the  angles  <J>1 ,  6^  = 

0,  6-j ,  62.  (b)  Specification  of  oxygen  locations  4  and  5  in  a  ring 

relative  to  the  dimer  via  bond  angles  d>2  and  plus  dihedral  angles 
(3124)  and  (2135).  Restrictions  on  6  for  9  >  120°. 

18.  A  3-oxygen  ring  with  the  silicons  in  the  plane  of  the  ring  illustrating 

I 

almost  total  constraint  of  the  6  's.  (b)  A  4-oxygen  ring  with 
0  =  160°  for  Si i  and  Si,,.  Here,  6  (Si ^ )  ~  0  to  40°,  and  6  (Si2) 

-  160°  ±  20°. 
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19.  Interaction  energy,  Ey,  of  02  with  the  dimer  as  a  function  of  sepa¬ 
ration  distance,  Z;  (a)  non-bonding  reaction  (02  in  plane  of  Si -0-Si ) , 
(b)  bonding  reaction  (02  perpendicular  to  plane  of  Si-O-Si). 

20.  Interaction  energy,  Ey,  of  H20  with  the  dimer  as  a  function  of 
separation  distance,  Z;  (a)  non-bonding  reaction  (H20  perpendicular 
to  the  plane  of  Si-O-Si),  (b)  bonding  reaction  (HO  in-plane  and  H 
perpendicular  to  the  plane  of  Si-O-Si). 

21.  Interaction  energy,  Ey,  of  H20  with  the  tetramer  for  a  number  of 
different  H20  configurations  in  the  ring.  Configuration  (g)  is  the 
most  favorable  (H-O-H  angle  =  105°  in  all  cases). 
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